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® Bar coding callbratioa 

0 For each test element being tested in a clinical analyzer, there Is a conresponding calibration curve. 
However, the calibration curve is defined by an equation having calibration coefficients which must be imparted 
to the analyzer so that the appropriate results can be calculated from the measurements taken off the test 
element. One way of imparting the calibration coefficients is to supply them in bar code fonm. However, due to 
the complexity of the coefficients, It may not be possible to store the information in a single bar code strip. 
Described herein Is a method of bar coding the data needed to determine a calibration curve for a test element 
in an analyzer by ascertaining the ranges of values for analyzer response (R) which are possible for three given 
concentration values (Ci . Cz, Ca), for a given assay, and assigning a high value (H) and a low value (L) for these 
ranges; caknjlating a calibration curve which correlates the analyzer response (R) to the concentration (C); 
determining the analyzer response (Ri , Rz, R3) from the calibration curve which conresponds to the concentra- 
tion values (Ci , O2, Cs); calculating the bar code value (B|) for the analyzer response (1%) where i = 1, 2, 3 from 
the equation: 

Bi = (10"-1)(R»-L4)/(H,-U) 

where U and H) are the corresponding low (L) and high (H) values for that Ri; rounding B| to the nearest integer; 
and supplying this value of B| for each of Ri in bar code form, using n-digit decimal numbers so that only three 
sets of (10") possibilities are needed. 
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This invention relates to the field of methods for calibrating a clinicaJ analyzer, and specifically to 
methods of passing on the data needed via a bat code to determine the calibration plot. 

It is conventional to calibrate a ctinical analyzer for a given assay and a given lot of test elements using 
several known calibration liquids with known analyte concentration (or activity). These liquids are dispensed 

5 on test elements from that lot, and responses are determined. The determined responses and the known 
concentrations are then used to compute calibration coefficients, using a known equation* so that such 
coefficients and equation can be used to calculate unknown concentrations using the responses generated 
from patient samples, using the same lot of test elements. 

For example, in a glucose test, it is conventional to use the equation 

10 Concentration = Ao + Ai.g^ (Response) + A2.g2(Response)'^ where g^ and g^ are cubic splines. K is an 
integer (usually = 2), and Ao, At and A2 are specific calibration coefficients. (See the "Principles of 
Caiibration'* section from the E700 Operators Manualj, The equation noted above has been published in 
connection with the analyzers available from Eastman Kodak Company under the registered trademark of 
"Ektachem". For simplicity, g(R) is hereinafter refen-ed to simply as "response", so that "R" is either the 

15 raw response or a cubic spline function of the raw response. 

It is also known that such calibration coefficients could be predetermined at the factory in some 
instances, and passed on to the purchaser of a given lot of test elements for that assay, to avoid making the 
user recalibrate each time a new lot is shipped. Such information is passed on in a variety of ways. 

The methods of passing on calibration information to the user include printed information and magnetic 

20 disks. Analyzers available from Eastman Kodak Company use a calibration diskette to transfer calibration 
infomiation to users but the calibration coefficients are not included on this diskette. It is not economically 
feasible to send a calibration diskette with each lot of slides for each assay, so calibration coefficients are 
not sent to the customer via a diskette. Lot specific canixation coefffcients can be transferred to the 
customer if the information is incorporated on the actual test element or its container. Two possible ways of 

25 transfening this information are by bar code or magnetic strip, as described in "Boelirlnger Mannheim 
detects high choiesteroi with the Reflotron diagnostics system". Directions, Vd. 6, No. 4. fourth 
quarter, 1989, or in Japanese Kokai 60/93351 . 

With single bar code strips and magnetic strips there is a limited number of digits, for example, six, 
available to pass the calibration coefficients, yet the purchaser requires accurate values for the calibration 

30 coefficients Ao. Ai and A2. If one is using a six decimal digit bar code, then Ao. Ai and A2 must be passed 
to the purchaser using only six decimal digits. The obvious solution is to specify Ac, Ai and Az using two 
digits each. This means that each of these must be accurately specified by using the digits 0 to 99, the 
maximum possible In a 2-digit decimal finite number. 

Prior to this invention, it has been a problem that the specific solutions of Ao. Ai and A2 can vary much 

35 more than this in a given assay, say glucose, because these coefficients are a function of the cubic splines 
that are used to best-fit the particular chemistry of a particular krt of test elements, to the data. Yet, bar 
coding is by far the preferred method of conveying the information of these coefficients, since that can be 
easily printed on each set of test elements or the package therefor. Nevertheless, it js well-recognized that 
a single strip of bar-coding is insufficient to portray the parameters of the calibration coefficients, as 

40 explained in Japanese Kokai 60^93351 . 

Of interest is the fact that the disclosure of Japanese Kokai 60/93351 attempts to solve the problem by 
providing, not a single strip of bar-coding, but rather, a triple strip of bar-coding, so as to allegedly Increase 
the number ol digits available to 1728 (12^). However, this approach is unsatisfactory since it requires both 
a much larger label for triple the amount of codes, as well as a much more sophisticated bar code reader. 

45 H is therefore an object of the present invention to provide a method of determining the data for the bar 
code which solves the problems mentioned above, while still using only a single bar code strip. 

More specifically, in accordance with one aspect of the present invention, there is provided a method of 
providing data in bar code form useful for the determination of the calibration curve of a lot of test elements 
in a clinical analyzer using a finite numbering system limited to n-digit decimal numbers, the curve having 

50 the mathematical fomoula:- 

C = Ao + Ai.R + A2.(R)'^ (1) 

where C is the predicted concentration of a sample liquid analyzed by the analyzer, R is the response 
56 actually measured in the analyzer or a cubic spline function of that response, K is a coefficient assigned to 
the analyzer, and Ao. Ai and A2 are the calibration coefficients which can vary beyond that which can be 
specified using (10") digits; 

characterized in that the method comprises the steps of:- 
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a) ascertaining by statistical analysis the ranges of values for R which are possible for three given 
concentration values (Ci. C2. Ca) for a given assay, and assigning a high value (H) and a low value (L) 
for these ranges: 

b) calculating for a given lot of the given assay, a calibration curve which correlates the analyzer 
5 response to the concentration; 

c) detenmin'mg the analyzer responses (Ri. R2, R3) from the calibration curve which correspond to the 
concentration values (Ci , C2. Ca); 

d) calculating the bar code value (B|) for the analyzer response (Ri) where i = 1, 2, 3 from the equation:- 
w B, = (10"-1)(R,.Li)/{H,.W (2) 

where U and H| are the corresponding low (L) and high (H) values for that value of R,; 

e) rounding the bar code value (B|) to the nearest integer for each case where \ = 1, 2, 3; and 

f) supplying this value of Bi for each of R| in bar code form so that only three sets of (10") possibilities 
f 5 are needed to accurately pass along data corresponding to the calibration coefficients even though each 

of the three coefficients can vary by more than that which can be specified using 10" digits. 
Accordingly, it is an advantageous feature of the invention that a single bar code strip of only a few 
digits can be accurately provided with the data needed to pass on a caTibration curve for a given lot of test 
elements, to the user. 

20 It is a related advantageous feature of the Invention that a simplified bar code, and hence a simplified 
bar code reader, can be used to represent the data needed to calibrate for lot-spedfic calibration 
parameters. 

For a belter understanding of the present invention, reference will now be made, by way of example 
only, to the accompanying drawings in which:- 
25 Figure 1 is a calibration plot of the expected concentrations for a given response in an analyzer for a 
given assay; 

Rgure 2 is an enlarged, fragmentary plot of just the abscissa axis of Rgure 1 ; 
Figures 3 and 4 are two-dimensional plots of a three-dimensional space of the possible values for the 
calibration coefficients Ao, At and A2 for the curves shown in Figure 1; and 
30 Figure 5 is a graphical representation of the differences which exist between the two methods shown in 
the Examples. 

The invention is hereinafter described in connection with the preferred embodiments, which use a 
preferred bar code on preferred dried, slide test elements in a preferred clinical analyzer. In addition, the 
invention is useful regardless of the form of the bar code, regardless of the format or assay of the test 
35 element (or Hs cartridge) on which the code is placed, and regardless of the analyzer in which the test 
element is tested. However, the invention is most useful in a single strip bar code. 

Any bar code design is useful with this Invention, provided that at least six digits are available, that is. 
three pairs, to provide three numbers which can range from 0 to 99. A particularly well-known form which 
provides this capability Is the so-called "interleaved two of 5". 
40 The preferred test elements are the slide test elements available from Eastman Kodak Company under 
the trademark "Ektachem" slides, or from Fuji Photo Rm Co. under the tradename "Drychem" slides. 

The prefeffed analyzers are any of the analyzers available from Eastman Kodak Company under the 
trademark "Ektachem" analyzer, or from Fuji Photo Rim Co. under the tradename "5000". 

Referring to Rgure 1, a representative plot is shown of a useful calibration curve, for example, for 
45 glucose. In such plots, the expected concentration C is plotted versus the response R measured on the 
analyzer, where R can be the raw response or a g(R) which is a cubic spline of the raw response. In 
general, the raw response can be any of reflectance, optical density obtained from reflectance, rate of 
change of these responses, or an electrical potential created by a differential measurement of ion 
concentration in two ion-selective electrodes. For glucose, the raw response Is either in reflectance or 
50 optical density Dr, where Dr = log(1 /reflectance). The curve can be expressed as concentration 

C = Ap + Ai .Response + A2 .Response*^ 

where K usually has a value of 2. 
55 The two plots, one a solid line and one a dashed line, both represent a good fit to the data that can be 
obtained on a given lot of test elements for this assay. That is. both curves occupy approximately the same 
space. However, the values of the coefficients Ao, Ai and A2 are drastically different for the two curves, as 
shown, where the response is optical density as determined by the analyzer. These values were determined 
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as follows: 

If one assumes the concentration C for three calibrators of different levels Is 39, 309, and 596mg/dl, as 
shown in Figure 1, and a corresponding response of 0.28, 1.1 and 1.5 Dr respectively, it is possible to solve 
for Ao, Ai arxJ A2 In the three linear equations ()), (II) and (lll)> 

5 

39 = Ao + Ai.gi(0.28) + P^.g2{02Qf (I) 

309 = Ao + A,.gi(1.1) + As.gzO.I)*^ (II) 

70 596 = Ao + Ai.g2(1.5) + A2.g2(1.5)'^ (III) 

A useful method for evaluating splines Is to use the spline parameters X, Y and F"(x), where P*(x) is the 
second derivative of the function at that x value, as described in "industrial Applications of Cubic Spllnei 
Functions", by N.J. Barosi, October 26, 1973, pp. 3-6 {A Presentation to tlie 17th Annual Technical 
IS Conference of the American Society for Quality Control and the American Statistical Association), 
and 'Splines and Statistics', by Edward J. Wegman and Ian W, Wright, Journal of the American 
Statistical Assdciation, June 1963, Volume 78, Number 382, Theory and Methods Section, pp. 351- 
352. 

20 Comparative Example 1 

Letting K = 2 and assuming the following spline parameters, where gi = g?:- 



25 



30 



X 


Y 


F"(x) 


-0.1 


0.0107 


0 


0.15 


0.2464 


0 


0.3 


0.4155 


6.54 


0.7 


1.422 


3.05 


1.2 


3.633 


6.65 


1.4 


5.012 


20.17 


2.0 


12.51 


0 



one finds that Ao, Ai and A2 are 19, 100.05 and 0.0025 respectively. See, for example, Rgure 1 for a curve 
representing these values (the solid line). 

Predicting the densities 0.2. 0.4, 0.6, 0.8, 1.2. 1.4, 1.6, and 1.8 using these calibration coefficierits, 
splines, and K-value one obtains concentrations of 30, 59, 111. 178, 260, 364, 502. 707. and 966 
respectively. 

Comparative Example 2 

K one changes the spline parameters so that they are 
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X 


Y 


F"(x) 


-0.1 


-0.0325 


0 


0.15 


-0.0166 


0 


0.3 


-0.0055 


0.436 


07 


0.0617 


0.203 


1.2 


0.2091 


0.444 


1.4 


0.3011 


1.35 


2.0 


0.8014 


0 



one finds that Ao, Ai and A2 are 50, 1500 and 0.5571 respectively. This is the dashed curve of Figure 1. 

Predicting the densrhes 0.2, 0.4, 0.6, 0.8. 1.0, 1.2. 1.4, 1.6, and 1.8 using these calitiration coefficients, 
splines, and K-value one obtains concentrations of 30, 59, 111, 178. 260. 364, 502. 707, and 966 
respectively. 

Thus the splines in Comparative Example 1 and Comparative Example 2 produce dramatically different 
calibration parameters Ao, Ai and A2. yet define the identical density-to-concentration relationship. Thus, 
coefficient A, for one curve is 1500, but for the other curve is 100.05. Clearly, any attempt to confine sucti 
variances of At to two digits of from 0 to 99 in value is doomed to failure, in terms of accuracy. Similar 
problems exist for bar^joding the variances obtainable just from these cun^es for Ao and A2. 

Still further, accuracy problems exist in bar coding the Ao. A, and Ag coefficrents. This Is illustrated In 
Comparative Example III which follows the description of the preferred embodiments. 

In accordance with the invention, instead of trying to fit the drastically varying Ao, Ai and Aa coefficients 
(of which Ao, Ai and Az noted above are specific examples) into the bar code, the solution is to fit the 
variances In the response R into the bar code. Such variances are in fact much less, as shown by the 
parentiiesis around Ra in Figures 1 and 2. The following non-exhaustive examples illustrate the practice of 
tiie invention. 

Example 1 (Glucose) 

The optical densities (Dr) associated with glucose concentrations of 40. 150 and 550mg/dl were found 
for five different generations of glucose slides. Data on twenty-one different coatings were found, using an 
"Ektachem 700" analyzer. The mean and standard deviation of tiie optical densities at the glucose 
concentrations were: 



Table I 



Concentration 


Mean Dr 


Standard Deviation (SD) 


40 


0.3122 


0.0594 


150 


0.7340 


0.0411 


550 


1.4765 


0.0381 



It can be shown, for a given concentration, that tfie optical density (Dr) on a new coating will fall, in 
99% of the cases, within the inten^al [Mean - 3SD. Mean + 3SD]. Thus, for the three fixed concentrations, 
new coatings must have optical densities for these concentrations which fall In the ranges: 
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Table II 



Concentration 


Low Dr Range (L) 


High Dr Range (H) 


40 


0.1399 


0.4905 


150 


0.6106 


0.8574 


550 


1.3623 


1.5908 



Therefore, for the equation (2) noted above, it is these values of H and L that are used to calculate what 
each B| should be for a given optical density (Dr) response at a given concentration C. 

More specifically, given a lot of slides, a glucose concentration C of 150mg/dl produces a Dr of 0.7969 
on one of the elements of the lot. It is this number which Is to be approximated using a two-digit barcode. 
Using the Dr ranges in Table II it is clear that the response of a fluid with a concentration of 150mg/dl must 
lie in the interval [0.6106. 0.8574]. Thus, the barcode value B2 is 

B2 = 99((0.7969 - 0.6106)/(0.8574 - 0.6106)).75 from equation (2). 

The approximate value found when the barcode is converted bacic by a customer using again an 
"EIrtachem 700" analyzer is: 

i^ (Converted) = 0.6106 + (75(0.8574 - 0.6106))/99 = 0.7976 

These converted values of Ru B2 and Ra are processed by the analyzer to create a new calibration 
curve using three sets of equations similar to equations (I), (II) and (III) above. In this case of R2. an error of 
only 0.0007 was induced by creating and converting the barcode. 

if one assumes that concentrations of 40 and 550mg/dl produce DrS of 0.4587 and 1.3706, this 
combination of concentrations and DrS. including a Dr of 0.7969 at 150mg/dl. produce a tme calibration 
curve with calibration coefficients of: 

-27.475, 68.737 and 6.764 using the spline of Comparative Example 1 above and K value of 2. 
Now, 

Bi = 99((0.4587 - 0.1399/(4.4905 - 0.1399)) 

= 90.02 X 90 (from equation (2)) 
B3 = 99((1.3706 - 1.3623)/(1.5908 - 1.3623)) 

= 3.60 X 4 (from equation (2)) 

and therefore, 

Ri (converted) « 0.1399 + ((90(0.4905 - 0.1399) )/99 
= 0.4586 

R3 (converted) = 1.3623 + ((4(1.5908 - 1.3623) )/99 
= 1.3715 



Calibrating with concentrations of 40, 150 and 550mg/dl and the converted densities of 0.4586. 0.7976. 
and 1.3715 one obtains the "converted" calibration coefficients: 
-27.227, 88.422 and 6.771. 

From these, a new curve is drawn and the densities corresponding to the "true" concentrations set forth 
in Table lit that follows, can be used to predict a "converted concentration". The difference between the 
"true" concentration and the "predicted" concentration using the "converted response" shows the error 
which would result from passing the above true calibration curve using the bar code which carries the 
response value. 
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Table III 





Tru6 Concentration 


Predicted Concentration UsinQ 


Absolute Bias 


5 




Converted Responses 






30 


30.06 


0.06 




61 


60.96 


-0.04 




92 


91.87 


-0.13 


10 


123 


122.79 


-0.21 




154 


153.70 


. -0.30 




185 


184.03 


-0.37 


15 


216 


215.55 


-0.45 




247 


246.48 


-0.52 




278 


277.41 


•0.59 


20 


309 


308.35 


-0.65 




340 


339.28 


-0.72 




371 


370.22 


-0.78 




402 


401.16 


-0.84 


25 


433 


432.11 


-0.89 




464 


463.05 


-0.95 




495 


494.00 


-1.00 


30 


526 


524.94 


-1.06 




557 


555.89 


-1.11 




588 


586.84 


-1.16 


35 


619 


617.80 


-1J20 




650 


648.75 


-1.25 



Such biases are negligible, as is seen in Figure 5. 

Although this example shows an "E700' analyzer being used both to create the bar code, and to adjust 

ttie calibration curve at the customer site, that need not be the case. That is, the site analyzer can be 
sirghtiy different from the one used to create the bar code, so long as the calibration mathematics is 
substantially the same for Iwlh types of analyzers. Thus, the bar code created on the "E700" analyzer can 
be used at a customer site which has an "E400" or "ESOO** analyzer, for example. 

^ Example 2: Blood Urea Nitrogen (BUN) 

The optical densities (Dr) associated with BUN concentrations of 10, 40 and 120mg/dl were found for 

two different generations of BUN slides. Data on sixteen different coatings were found using an "Ektachem 

E700". The mean and standard deviation of the opticai densities at the BUN concentrations were: 
so ^ 
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Table IV 



Concentration 


Mean Dr 


Standard Deviation (SD) 


10 


0.3608 


0.0231 


40 


0.7958 


0.0278 


120 


1.8725 


0.0755 



It can be shown, for a given concentralion, that Dr on a new coating will fall within the interval [Mean - 
3SD, Mean + 3SDJ. Thus, for the three fixed concentrations, new coatings must have DrS for these 
concentrations which fall in the ranges: 



Concentration 


Low Dr Range (L> 


HIgti Dr. Range (H) 


10 


0.2915 


0.4301 


40 


0.7124 


0.8792 


120 


1.6458 


2.0991 



Therefore, for the equation (2) noted above, it Is these values of H and L that are used to calculate what 
25 each B| should be for a given Dr response at a given concentration C. 

More specifically, given a lot of slides, a BUN concentration of C of lOmg/dl conresponds to a Dr of 
0.3756 on one of the elements of the lot It is this number that is to be approximated usin9>a two digit bar 
code. 

Using the Dr ranges In Table V it Is clear that the response of a fluid with a concentration of lOmg/dl 
30 must lie In the interval [0.2915, 0.4301]. Thus the bar code value Bi Is 
Bi = 99 ((0.3756 - 02915)/(0.4301 - 0291 5))60 from equation (2). 

The approximate value found when the bar code is converted back by a customer using again an 
"Ektachem 700" analyzer Is: 

Ri (converted) »0.2915 + ((60(0.4301-0.291 5)y99K0.3755) 
35 In this case, an error of only 0.0001 was induced by creating and converting the bar code (an error of 
only 0.03%). 

ft will be appreciated from the foregoing examples that the bar-coding of the particular response Ri as a 
fraction of (H - L) for given values of Ci, and Ca, acts to bypass a passing on of the coefficients Ao, Ai 
and A2. Ao, Ai and A2 are instead recomputed in the analyzer using the Fli and Ci as described above. 

40 

Comparative Example ilk Barcoding Calibration Parameters Ao, Ai and A2 

Figures 3 and 4 indicate another weakness in specifying calibration parameters rather than test element 
responses. In this study the following simulation was performed: 
45 Step 1) For each of the glucose concentrations in Table II, that Is, 40, 150 and 550, the range listed was 
partitioned Into ten evenly spaced densities which spanned the known density ranges for the concentration. 
The resulting densities for the given concentrations are listed in the Table VI below. 



50 
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Table VI 



10 



15 



40mg/dl 


ISOmg/dl 


550mg/dl 


0.1399 


0.6106 


1.3623 


0.1789 


0.6380 


1.3877 


0.2187 


0.6854 


1.4131 


0.2568 


0.6929 


1.4385 


0.2957 


0.7203 


1.4639 


0.3347 


0.7477 


1.4892 


0.3736 


0,7751 


1,5146 


0.4126 


0.8026 


1.5400 


0.4515 


0.8300 


1.5654 


0.4905 


0.8574 


1.5908 



20 

Step 2) K was set equal to 2 and gi = g2 = glucose spline from Comparative Example 1 at>ove. 

Step 3) All possible combinations of densities from Table V were found where one density was chosen 
from each column. This resulted In one thousand sets of three densities. 

Step 4) One thousand calibrations were performed using the glucose concentrations 40, 150 and 
25 550mg/dl and each set of densities found In step 3. 

Step 5) The one thousand sets of calibration parameters Ao. Ai and A2 generated in Step 4 were 
recorded. 

Step 6) Ao vs A? from each calibration parameter set found in Step 5. was plotted in Figure 3. and Ai 
vs A2, from each calibration parameter set found in Step 5, was plotted in Rgure 4. (Plotting all three 
30 parameters at the same time would require creating a three dimensional ploL) 

Now, if one attempted to bar code Ao, Ai and Az rather than the test element responses, one would 
need to specify three ranges using the same format as In Table II for Aq, Ai and A2. These ranges would 
be similar to the ones listed in Table Vli below, detenr^ined from Figures 3 and 4: 

35 Table Vli 



40 



Ao» Ai and A2 


Low Range 


High Range 


Ao 


•300 


100 


Ai 


0 


500 


Az 


-60 


20 



When bar coding the True Calibration Coefficients (that is. -27.475, 88.737, and 6.764) found above 
using Table VII, one obtains: Bi = 99((-27.475 + 300)/(100 + 300)) = 67.45 x 67 B2 = 99((88.737 - 0) / 
(500 - 0)) = 17.57 X 18 B3 = 99((6.764 + 60) / (20 + 60)) = 82.62 x 83 

Converting the bar code bacl( to calibration coefficients one finds: 
Aq (converted) = -300 + (67(100 + 300))/99 = -29.293 
Ai (converted) = 0 + (18(500 - 0))/99 = 90.909 
A2 (converted) = -60 + (83(20 = 60))/99 = 7.071. 
Thus, the Converted Calibration Coefficients are -29,293. 90.909. and 7.071. 

Table Vill below shows the error in predicted concentration which would result from passing the true 
calibration curve using calibration coefficients on the bar code, rather than the responses as per the 
invention. 
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Table VIII 





ilue \>onceniiaiion 


riuuicicGi Vroncsnuraiion usinQ 


ADsoiuxe Dias 


5 




Converted Responses 






30 


29.64 


-0.36 




61 


61.47 


0.47 


70 


92 


93.32 


1.32 


123 


125.19 


2.19 




154 


157.08 


3.08 




185 


188.99 


3.99 


16 


216 


220.91 


4.91 




247 


252.84 


5.84 




278 


284.79 


6.79 


20 


309 


316.75 


7.75 




340 


348.72 


8.72 




371 


380.70 


9.70 


25 


402 


412.70 


10.70 


433 


444.69 


11.69 




464 


476.70 


12.70 




495 


508.72 


13.72 


$0 


526 


540.74 


14.74 




557 


572.77 


15.77 




588 


604.81 


16.81 


35 


619 


636.85 


17.85 
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Figure 5 Is a graphical representation of the differences which exist between the two methods for this 
^ example. It uses the data found in tables III and VIII. (The ordinate value of zero represents "truth".) A 
comparison of Table III and Table Vltl shows that passing responses on the bar code results In a 
significantly better approximation of the true calibration curve than the calibration curve found by passing 
calibration coefficients directly on the bar code. 

^ Claims 

1. A method of providing data in bar code form useful for the determination of the calibration cun^e of a lot 
of test elements in a clinical analyzer using a finite numbering system limited to n-digit decimal 
numbers^ the curve having the mathematical formula:- 

C = Ao + Ai.R + A2.(Rf (1) 

where C is the predicted concentration of a sample liquid analyzed by the analyzer. R is the response 
actually measured in the analyzer or a cubic spline function of that response, K is a coefficient 
^ assigned to the analyzer, and Ao. Ai and A2 are the calibration coefficients which can vary beyond that 
which can be specified using (10") digits; 

characterized in that the method comprises the steps of:- 

a) ascertaining by statistical analysis the ranges of values for R which are possible for three given 
concentration values (Ci, C2. Ca) for a given assay, and assigning a high value (H) and a low value 
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(L) for these ranges; 

b) calculating for a given lot of the given assay, a calibration curve wtiich correlates the analyzer 
response to the concentratiori; 

c) determining the analyzer responses (Ri. R2. R3) from the calibration curve which con'espond to 
the concentration values (Ci , Ca. C3); 

d) calculating the bar code value (Bi) for the analyzer response (R|) where i - 1. 2, 3 from the 
equation:- 

Bi = (10"-1)(R|-l4y(H,-l^) (2) 

where U and Hj are the corresponding low (L) and high (H) values for that value of Rj; 

e) rounding the bar code value (Bi) to the nearest integer for each case where i - 1, 2, 3; and 

f) supplying this value of B| for each of R| in bar code form so that only three sets of (10") 
possibilities are needed to accurately pass along data corresponding to the calibration coefficients 
even though each of the three coefficients can vary by more than that which can be specified using 
10" digits. 

A method according to claim 1, further including the step of converting the supplied bar code values 
(B|) of step e) into actual values of analyzer responses (R|) by solving equation (1) for Ri, given the 
stored values of B|. 

A method of providing data according to claim 1 or 2, wherein n = 2, and the step a) includes the 
steps of determining the mean value (M) to be expected in all test elements for the given assay and the 
standard deviation (SD) from that mean (M), and assigning the high value (H) and the low value (L) as 
M + 3SD and M - 3SD respectively, for a given concentration value (Ci » Ca, C3). 
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